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Abstract 

We study low-energy properties of the random displacement model, a random Schrddinger 
operator describing an electron in a randomly deformed lattice. All periodic displacement con- 
figurations which minimize the bottom of the spectrum are characterized. While this configu- 
ration is essentially unique for dimension greater than one, there are infinitely many different 
minimizing conhgurations in the one-dimensional case. The latter leads to unusual low energy 
asymptotics for the integrated density of states of the one- dimensional random displacement 
model. For symmetric Bernoulli-distributed displacements it has a 1/ log 2 -singularity at the 
bottom of the spectrum. In particular, it is not Holder-continuous. 

I. Introduction 

We consider the so-called random displacement model, i.e. the random Schrddinger operator 

H u = —A + V u , (1) 

where the random potential V w is given by displacing a single site potential q from the points 
of 7J ! . 

V u (x) = ^2q(x-i-u}i). (2) 

For the real- valued single site potential q we assume q G and suppg C [— r, r] d for 

some r < 1/2. We also assume that q is reflection symmetric at each coordinate hyperplane, 
i.e. symmetric in each variable with the remaining variables fixed. Throughout this paper we 
will consider displacement configurations u = (uJi) ieZ d such that Ui G [—d max ,d max ] d for all i, 
where r + d max = 1/2. The latter ensures that the displaced single site potentials in (j2J) do 
not overlap. 

While the random displacement model is a physically quite natural way to describe struc- 
tural disorder in a solid, it is mathematically much less well understood than Anderson-type 
models, where the disorder enters in the form of random couplings at the single-site potentials. 
This is mostly due to the fact that the random displacement model depends non-monotonously 
on the random parameters (in quadratic form sense). Anderson- type models, on the other 
hand, are monotonous in this sense, at least if the single-site potentials have fixed sign. The 
consequential challenge in determining the spectral properties of the random displacement 
model, in particular the low energy behavior, lies in having to gain a deeper understanding of 
the interaction between multiple random parameters as well as the interplay between kinetic 
and potential energy. Spectral averaging arguments using individual random parameters, a 
common tool in the theory of Anderson models, are not available here. 
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This is one of the main reasons why it is not yet known if the multi-dimensional random 
displacement model is localized at the bottom of the spectrum (in d — 1 localization at all 
energies follows from the results in [6], see also [T8]). The only known result on localization for 
the multi-dimensional random displacement model is due to Klopp [TO], who considered the 
semi-classical version —h 2 A + V w of the random displacement model and identified a localized 
region near the bottom of the spectrum for sufficiently small values of the semi-classical 
parameter h. 

An attempt to understand the low energy properties of the random displacement model has 
to start with describing the mechanism which characterizes the bottom of the spectrum. Under 
the above assumptions this has been achieved in [2] by identifying the periodic displacement 
configuration u min , see 01]) below, which leads to the minimum of the almost sure spectrum 
of Hu. 

Here we continue our study of the low energy properties of the random displacement model 
by first characterizing the set of all minimizing periodic configurations. It turns out that u) mm , 
up to trivial translations, is the unique minimizer in d > 2, while in d — 1 there are infinitely 
many periodic minimizers, which can be explicitly characterized. These results are stated in 
Section [2] and proven in Section [3j 

We then move to studying the low energy asymptotics of the integrated density of states 
(IDS) for the random displacement model. Showing smallness of the IDS near the bottom of 
the spectrum, typically in the form of Lifshits tails, is an important step in all approaches to 
low-energy localization for multi-dimensional random Schrodinger operators. It is interpreted 
as showing that the bottom of the spectrum is a fluctuation boundary. 

However, as we will show here, for the one-dimensional random displacement model the 
behavior of the IDS can be very different. If the displacements only take values d max or 
—dmax, both with equal probability, then the IDS N(E) has a very strong singularity at the 
bottom of the spectrum, 

log (E-Eq) 

for E e (E Q , E Q + e) and constants C > and e > 0, see Theorem 14.11 Here E denotes the al- 
most sure minimum of the spectrum of H u . Thus the IDS is not even Holder-continuous at Eq, 
a new phenomenon which to our knowledge has not been found for any other models of random 
operators. This and related results for the IDS of the one-dimensional random displacement 
model with other distributions of the displacements are proven in Section HI For example, the 
extreme behavior ([3]) only appears for the given case of a symmetric Bernoulli distribution 
of the displacements, see Theorem 14.21 But, in d = 1 and as long as the distribution of the 
displacements u n is symmetric, one never gets Lifshits tails (Theorem 14. 3p . 

In Section [S] we comment on several open problems, in particular on our expectation that 
the uniqueness of the periodic minimizing configuration in d > 2 should indicate different low 
energy asymptotics for the IDS than in d = 1 (e.g. the appearance of some form of Lifshits 
tails). We will also discuss related recent works by Klopp and Nakamura [11] and Fukushima 

2. Periodic configurations which minimize the ground state energy 

In [2] we have identified a simple periodic configuration of displacements which leads to the 
lowest possible spectral minimum for H w among all configurations u: 
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Proposition 2.1 (Theorem 1.1 in |2j). Let u mm be given by 

^min = ((_l)*i dmaxj . . . , (-iy*d max ) (4) 

for all i = (ii, . . . ,id) G Z d . Then 

E Q = mincr(if w min), (5) 

where E := inf^mincr(i^ w ). 

This configuration is 2-periodic in each coordinate, where in each period cell 2 d single sites 
cluster together in adjacent corners of unit cubes, see Figured] for d = 2. 




Figure 1. The support of V w min for d = 2. 

Here our first goal is to characterize all periodic configurations uo such that min ^{H^) = Eq. 
For this we will also use another result found in [2]: Let q be as above, a e [—d max ,d max ] d , 
A = (— J, \) d the unit cube centered at 0, and H^ a (a) = —A + q(x — a) on L 2 (A ) with 
Neumann boundary condition on OAq. 

Proposition 2.2 (Theorem 1.3 of [2J). // 

E (a):=mma(Hl(a)), (6) 

then the following alternative holds: Either 

(i) E (a) is strictly maximized at a = and strictly minimized in the 2 d corners 

maxi • • • j i^max) of [ d max , d max ] , Or 

(ii) E (a) is identically zero. In this case the corresponding eigenfunction is constant 
outside of the support of q. 

Alternative (i) holds for generic potentials, for example if q is non-zero and sign-definite, in 
which case Eq(o) never vanishes. But sign-definiteness of q is far from necessary for alternative 
(i). Alternative (ii) holds if and only if the Neumann problem for — A + q restricted to (— r, r) d 
has ground state energy (if suppg is simply connected with sufficiently regular boundary 
one may equivalently use the Neumann problem on supp q here) . Examples of sign- indefinite 
q with this property can easily be constructed. 

Theorem 2.3. (a) If alternative (ii) holds, then mincr(if aj ) = for all configurations u = 
{Ui) i& d, with uji G [-d max , d max ] d for all i. 

(b) If alternative (i) holds, d>2 andr < 1/4, then uj mm given by O) is, up to translations, 
the unique periodic configuration with min a{H UJ min) = Eq. 
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We don't believe that the extra condition r < 1/4 (beyond r < 1/2) is required in Theo- 
rem 12.31 but we need it in our proof. 

It remains to settle the case of alternative (i) and d — 1. In this case the periodic minimizer 
is highly non-unique, but the set of all periodic minimizers can be characterized by our next 
result. Here, for L G N let Sl denote the set of all L-periodic configurations (ui) ie i such 
that Ui = —d max or uji = d max for all i. Furthermore, for u G Sl let n (u) be the number of 
z G {1, . . . , L} with Ui = ±d max . 

Theorem 2.4. Let d = 1 and q such that alternative (i) holds. An L-periodic configuration 
u satisfies min o~(Hj) = E if and only if L is even, u G Sl, and n~(u) = n + (u). 

Thus, in each period interval of V w , equally many of the single site potentials sit at the 
extreme right and the extreme left of their allowed range of positions. The dimer configuration 
u} mm is merely a special case of this situation. 

The short explanation behind Theorems 12.31 and 12.41 is that for d = 1 the ground states of 
H Ao (a) for extremal positions of a can be made to match after suitable re-scaling, while the 
richer geometry for d > 2 prevents this for all configurations other than u mm . 



3. Proof of Theorems 12.31 and 12.41 

Part (a) of Theorem 12.31 easily follows from known facts: Under alternative (ii) the ground 
state energy of the Neumann problem for —A + q restricted to (— r, r) d is 0. Let <p be the 
corresponding positive normalized eigenfunction. For any given configuration u place a trans- 
late of (p at i + out for each i G Z d and extend by a constant to the exterior of the convex 
hull of suppg(- — i — Ui). This results in a globally bounded, positive weak solution (p of 
— Aip + V^(p = 0. By Shnol's Theorem we have G a(H u ). As (p is positive it also follows 
from Theorem C.8.1 in [14J that < min ^(iL^). Thus wxx\.<j(H<J) = for all configurations 
u. 

For the remainder of this section we will assume that alternative (i) holds. If u is a periodic 
configuration with period L = (Li, . . . , La) G N d , i.e. such that Wt+(niii,... ) n d L <J ) — f° r ah i £ 
Z d and and (m,...,n d ) G Z d , we choose A = (1/2, . . . , Li + 1/2) x . . . (1/2, . . . , L d +l/2) C R d 
as period cell of the potential V u given by (T5]). By H^ A and H^ A we denote the restrictions of 
— A + K, to L 2 (A) with periodic and Neumann boundary conditions, respectively, and denote 
their lowest eigenvalues by E (H^ A ) and E (H^ A ). 

The proofs of Theorems 12.3( b) and Theorem 12.41 will be based on the following result 
which holds for arbitrary dimension. It shows, in particular, that in minimizing periodic 
configurations all single site potentials necessarily must sit in the corners of unit cubes centered 
at the points of Z d . 

Lemma 3.1. Let u be a periodic configuration with mino^if^) = Eq. Then, for all i G 
Ij d , Ui G {(ai, • • • , a>d) G M d : afc G {—d max , d max } for all k = 1, . . . , d}. Moreover, in this 
case E (H^ A ) = E (H^ A ) and the ground state eigenfunction ip^ of A satisfies Neumann 
boundary conditions on the boundary of each unit cube Aj centered at i G A D 

Proof. By assumption and Floquet-Bloch theory Eq = mincx(if a; ) = E (H^ A ). Also, by the 
variational principle, E (H^ A ) > E (H^ A ). The ground state ip u minimizes the quadratic 
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form of H^ A , thus 



/ A |V^| 2 + f A EieAnz^ Wi)|^ a 



E 



J Aj + Xv, ^)hU 2 J A . \A 



> 



E^)f^> E^of^^o, (T) 

where -Eo( w i) is given by (jSJ)- In the second to last inequality we have used the variational 
principle as well as the fact that if>^ does not vanish on any of the Aj. 

We conclude that all inequalities above must indeed be equalities, which immediately gives 
E o{H^ A ) = E (H^ A ). If for at least one Ui G" {a : a k G {-d max , d max }, k = 1, . . . , d}, then the 
last inequality in would be strict, given that we are in alternative (i). We conclude that 
all Ui sit in a corner. Finally, we see that 

12 

- = E (ui) (8) 



/ A . \V^\ 2 + j^x-i-u^ 



for each i. Thus the restriction of if> u to Aj is the ground state for the Neumann problem of 
— A + V(x — i — Ui) on Aj and thus satisfies Neumann boundary conditions on each Aj. □ 

We now consider d = 1, still under alternative (i), where we use one more lemma to prepare 
for the proof of Theorem 12.41 

Lemma 3.2. Let d = 1 and H Ao (a) be the restriction of —d 2 / dx 2 + q(x — a) to L 2 (— |, |) with 
Neumann boundary conditions. Let if) be the positive normalized ground state of H Ao (d max ). 
Thenif){±)^if){-\). 

Proof. Suppose that ip(\) = 4>(— \), then if) coincides with the periodic ground state of 
—d 2 /dx 2 + q(x — d max ) on L 2 (— |, |). Due to the symmetry of q and the uniqueness of the 
periodic ground state, it must therefore also satisfy if)'(—\ + d max ) = (considering (— |, |) 
as a 1-torus, — \ + d max lies opposite to d max ). As q(x) = for x G (— |, — | + d max ), we have 
that -if)" = E (d max )if> on [—\,-\ + d max ] and 

^{-\) = ip'{-\ + d max ) = V. (9) 

If E (d max ) < 0, then if)" = —E (d max )ip > and thus if) is strictly convex on [—1/2, —1/2 + 
dmax], contradicting (jUj). Similarly, Eo(d m ax) > would yield strict concavity of if), again 
contradicting ([9]). Thus E (d max ) = and it follows from (jSJ) that if) must be constant outside 
the support of q(x — d max ). This contradicts that we have assumed alternative (i). □ 

We can now complete the proof of Theorem 12.41 Let u be an L-periodic configuration which 
satisfies mino^if^) = E . By Lemma [3.11 we have u ri G {±d max } for all i = 1, . . . ,L and also 
E = rnin<7(#£ A ) = mina(H^ A ), where A = (|,L+ §). 

Let «o and w^v be the solutions of — u" + Vu = E Q u which satisfy = 0, u' D {\) = 1, 

un{\) = 1) m at(|) — 0. Eq is a Neumann eigenvalue on A and thus u' N (L + |) = 0. The 
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transfer matrix for at E from 1/2 to L + 1/2 is given by 

u N (L + l) u D (L+±] 



T -{u' N (L + l) u' D (L+l)J- ^ 
This implies 

1 = detT = u N {L + l/2)u' D (L + 1/2). (11) 
Moreover, as Eq is also an eigenvalue for periodic boundary conditions, 

2 = trT = u N (L + 1/2) + u' D {L + 1/2). (12) 
We conclude from ( flT|) and ( jl~2l) that 

Ujv (L + 1/2) = 1 = i4(1/2), (13) 

meaning that is both the Neumann and periodic eigenfunction to E . 

We can use Lemma [3721 to understand the detailed structure of mat: Let if) be the normalized 
ground state of H¥(d max ) as given there. Then, by symmetry of q, the normalized ground 
state of H^ Q (-d max ) is given by $(x) = ip(-x). As w» G {±<i max } for alii = 1, . . . , L, we 
can construct by concatenating suitably re-scaled versions of ip and ip, respectively, on the 
intervals [i — |,i + |]. With the positive number r = ip(^) |) 7^ 1 from Lemma [3.21 we 
thus have 

r if CO,- = d r 



u N {i+i/2)/u N {i-i/2) = \ \ -Tr ( 14 ) 

^ r 11 '-''max- 

The accumulative effect is that u N (L + §) = W^""-^^!). We conclude from ([T3J) that 
n + (to>) = n_(c<j), in particular that L is even, which completes the proof of Theorem 12.41 

We now start with preparations for the proof of Theorem 12.3( b). 

Lemma 3.3. Consider a connected open region D in IR d ; d > 2 and a hyperplane P that 
divides this region into two nonempty subregions. Denote by a the reflection about P and 
assume that D n er(-D) is connected. Let E e R and, in D, let u be a solution of the equation 

-Au = Eu (15) 

which satisfies the condition |^ = on PDD. Thenu can be extended to a symmetric function 
w on D U er(-D) which satisfies the equation —Au = Eu in this region. 

Proof. Pick a point xq G P D, which we may assume to be the origin. Pick a ball B C 
D centered at the origin and pick coordinates xi,...,x n so that x' = (x±, . . . , x n _i) are 
coordinates in P and x n is the coordinate normal to P. Consider the function 

which satisfies (1151) and vanishes on B D P identically. Its first normal derivative satisfies 

<9t> 



<9x, 

Since 



-(x',0) = 0. (16) 



= Evix', 0) = -Av(x', 0) = 0) (17) 

we also have that 

g(x',0) = 0. (18) 
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Further, since 







E 



( 




n J V C72'n / 



(19) 



(20) 




(x',0) = . 



(21) 



dxl 

Continuing in this fashion we deduce that 



(x',0) = 0, 



(22) 




for A: = 0,1,2,.... In particular, this implies that all derivatives of v vanish at the origin. 
Since v solves (TTo^) it is a real analytic function and hence it vanishes in the ball B. Since 
Xq G D fl <j(D) and D (1 &{D) is connected we learn that v vanishes everywhere in D D &(D). 
Thus u is symmetric in .D fl (?(D) with respect to reflection about the plane P. Next we 
prolong u to the complement of D fl cr(-D) in DU cr(-D) by setting 



Note that this function is defined since u(x) and u(a(x)) coincide on D fl cr(-D)- Moreover, 
x G D U cr(D) means that s G -D or x G er(.D) or both. In any case, by assumption and the 
fact that the Laplace operator commutes with reflections, w(x) satisfies the equation ffTBT) at 



Given the previous lemma, we can now complete the proof of Theorem 12.3( b). Suppose 
that d > 2 and u is a periodic configuration with mm.a{H UJ ) = E , but not a translate of 
uj mtn . Then there must be two adjacent unit cubes, say A and A', such that the potential on 
the union R = A U A' of these cubes is not symmetric with respect to reflection about their 
common face. This common face defines a hyperplane P. By Lemma 13.1} the two potential 
sites of Vu restricted to R are supported in corners of A and A', respectively, in the sense that 
[ai — r, at + r] x . . . x [ad — r, + r] sits in a corner of A, and similar for a' and A'. Also, the 
positive ground state eigenfunction ip^ of satisfies Neumann conditions on OR as well as 
on P. Let 

D = R \ ([ai — r, a x + r] x . . . x [a d — r, a d + r] U [a[ — r, a[ + r] x . . . x [a' d — r, a' d + r]) (24) 

and u the restriction of ip^ to D. For these choices of D and u and E = Eq we can apply 
Lemma 13.31 In particular, the assumption r < 1/4 assures the required connectedness of D 
and D fl cr(D). Non-symmetry of a and a' and again r < 1/4 implies that D U cr(D) is all of 
R. Thus, by Lemma 13. 3^ u can be extended to a symmetric function w on R which satisfies 
—Aw = Eqw on R. Therefore it is the ground state for the Neumann problem of —A on R. 
This implies E = and that w is constant, a contradiction to the assumption of alternative 




x G D, 
x G cr(D). 



(23) 



this point which proves the lemma. 



□ 
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4. Consequences for the integrated density of states 

We now consider the random displacement model, i.e. the model (OQ), ([2]) for the case where 
uj = (co>j)j gZ d is an array of i.i.d. vector-valued random variables with common distribution /x 
supported on [— d max , d max ] d . Here, as usual, we define 

supp fi := {a E R d : /x({x : \x — a\ < e}) > for all s > 0}. 

Then the random operator H w is ergodic with respect to translations in Z d and thus has all 
the basic properties of ergodic operators, see e.g. [3j. In particular, the integrated density of 
states (IDS) 

N(E) = lim -^-E(tT X (-oo,E](H^ AL )) (25) 

exists for all energies Eel. Here A L = (|,L + and H* Al is the restriction of H u to 
L 2 (Ax,) with boundary condition X 6 {P, N, D}, as periodic (P), Neumann (N) and Dirichlet 
(D) boundary conditions all give the same limit in ( 1251) . 

The spectrum a(H w ) is almost surely deterministic, i.e. E = o-(Hj) for almost every u, and 
given by the growth points of the non- decreasing function N(E). It can be characterized in 
terms of the spectra of those H w for which the configuration u is periodic, 

E = \Ja(H u ), (26) 

U! 

where the union is taken over all periodic u such that Ui e supp /x for all %. This corresponds 
to a well-known result for the Anderson model, e.g. [3], and is found with the same proof. If 
the support of the distribution /x contains all the corners {(±d max , . . . , ±d max )} of the cube 
[—d m ax,d m ax\ d , then it follows from (1261) and Proposition 12. II that 

minS = Eq = min a{H UJ m.i n ). 

For large classes of random Schrodinger operators it is known that the IDS vanishes rapidly 
at the bottom of the spectrum Eq, for example one has Lifshits tail behavior 

N(E) - e -\ E -^\- m (27) 

for Anderson models with sign-definite single-site potential, see e.g. [9], [13j0-ppp or [T5] for 
proper statements, proofs and references to the original literature. 

It turns out that for the random displacement model the behavior of the IDS at the bot- 
tom of the spectrum is much more subtle. Here we will present several results for the one- 
dimensional displacement model, which were obtained in pQ. We will generally assume that 
alternative (i) holds. 

In the first result we will consider the one-dimensional Bernoulli displacement model, i.e. 
the case where the distribution of the displacements u>j is given by 

1 1 

u = -5 d + -cL d . (28) 

It turns out that in this case the low-energy asymptotics of the IDS is at the opposite extreme 
of Lifshits tails: 

Theorem 4.1. Let be the one- dimensional symmetric Bernoulli displacement model given 
by & an d WW an d assume that alternative (i) holds. Then there exist C > and e > 
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such that 

N(E) > ° (29) 

log (E-Eq) 

forEe(E ,E + e). 

As N(Eq) = 0, this means that N(E) has infinite upper derivative at E — Eq, i.e. the 
density of states n(E) = N'(E) has a strong singularity at the bottom of the spectrum. This 
is opposed to the case of Lifshits tails which would yield n(E ) = 0. In fact, fl29|) says that 
the IDS is not even Holder continuous at E — Eq, an even stronger singularity than one gets 
for the Laplacian Hq = —d 2 /dx 2 , where the IDS has a van Hove sing ularity C\E\ 1 / 2 . For 
general one-dimensional ergodic Schrodinger operators (and for discrete ergodic Schrodinger 
operators also in higher dimension) the IDS is log-H61der-continuous at all energies, i.e. 

W(m-N(B)\< ^-B\\ [W) 

for E close to E', see (HE]. Craig and Simon constructed examples of quasi-periodic potentials 
which show that the bound (I30p is optimal. As far as we know, the result in Theorem 14.11 
provides the first known example of a random potential (with finite correlation length) where 
for at least one energy the IDS is not Holder-continuous and, in fact, close to the minimal 
possible regularity for ergodic operators given by (1301) . 

Proof. We will use the standard lower bound, e.g. [3], 

N(E) > i-P(£ « L ) < E), (31) 

which holds for arbitrary L, to be chosen later depending on E. Here H D L is short for H® A , 
A L = (l/2,L+l/2). 

To show that E (H D L ) < E we will find ^ e D(H° L ) with \\if> u \\ = 1 and (ip u , H^J) < 
E. To construct ip^, let displacements uj = (ui, . . . ,ul) be given and let ujy be the solution 
of —u" + V w u = E u with Ujv(|) = 1, u' N (\) = 0. Choose cut-off functions 6 L e C™(R) 
with <6 L <1, supp6 L C [1,L], 6 L (x) = 1 for 3/2 < x < L - 1/2, and H^H^ and \\6'l\\oo 
uniformly bounded in L. 

As the uji have distribution fT2S|) . we have uoi e {±d max } for all i. Thus the restriction of 
—d 2 /dx 2 + Vu to (i — 1/2, i + 1/2) has Neumann ground state energy E for all i G {1, . . . , L}, 
which implies that 

u' N {i + 1/2) = for aU % e {1, . . . , L}. (32) 
We choose tp u :— 9lUn/\\9lUn\\ and calculate 

(9 L u N ,-9lu N )-2(e L u N ,e' L u' N ) 



(^H^M-Eo 



< 



(3(1 +u%(L + 1/2)) 
I3/2 u n( x ) dx 



K P(l+u 2 N (L + l/2)) 

" £ti«3U< + i/2) 



where (3 > and j3 > can be chosen uniformly in uj and L. Here we have repeatedly used 
standard a priori upper and lower bounds on solutions of — u" + Vu = Eu, for example that 
r(x) ~ r(x + 1) and J x+1 u 2 ~ r 2 (x), where r(x) = (u 2 (x) + u' 2 (x)) 1 / 2 is the Priifer amplitude 
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of u and constants can be chosen uniform as long as E and ||V||oc vary in a bounded interval, 
see e.g. [HI [T7] for more details. We also use that by (J3"2"j) the Priifer amplitude of at the 
points 2 + 1/2 coincides with un (i + 1/2). 
Thus 

P(£o«J < £) > P + < E ~ A- (34) 

Another consequence of cjj G {±c/ ma:r } is that un satisfies ( JT4l) for every z with a positive 
r / 1, using that we are in alternative (i). Assume without restriction that r>l(ifO<r<l 
then we can do the following construction from "right to left", choosing un{L + 1/2) = 1, 
u' N (L + 1/2) = 0) and set 

Xi = log(M* + l/2)/u N (i - 1/2)) 
log r 

i — 1, . . . , L. The Xi are independent symmetric Bernoulli random variables with values ±1, 
and 

u 2 N {i + l/2) = e 2S ^% (36) 

where Si = X 1 + . . . + X^. If Y := max i=lv . jL Si, then it is a consequence of the reflection 
principle for symmetric random walks, e.g. [7] that 

¥(Y > VI\S L < 0) = F(S L > 2VZ). (37) 

The latter converges to 7r -1 / 2 f 2 °° exp(— y 2 /2) rfy > as L — > oo by the central limit theorem. 

Let A L := {lu\Y > Viand S L < 0}. If Y > y/I, then J^Li u n( l + V 2 ) > exp(2 v / Ilogr). 
Also, S L > means u 2 N (L + 1/2) < 1. Thus (JMD implies 

P( J E « L )<E) > P 1 ^ ^ + yfol <^-^o|^ ) HA L ) 

E t =i W + 1/2) 



1l) > c > (38) 

if 2/3 exp(— 2 v^log r) < E — Eq and L sufficiently large. This determines the choice of L G N 
for given E such that 

J_ -VE^liogr > E-E > ^ e - 2VZlogr . (39) 
2/3 - 2/3 v 7 

Thus L ~ ( 1 ° g2 ^~ i;o) ) 2 - From (EU) and ([38]) we have iV(£) > cq/L, which, for E - E 
sufficiently small, takes the form (1291) . 

□ 

As mentioned above, (1291) says in particular that the IDS is not Holder continuous at 
E = Eq. This is only possible if the distribution \x is concentrated in the extreme points d max 
and — d max , as is demonstrated by our next result. 

Theorem 4.2. Suppose that the distribution fi of the u>i in the one-dimensional displacement 
model (TJP, (0|) satisfies 

M(- dmax, d max )) > 0. (40) 
Then the IDS N(E) is Holder continuous at E = Eq. 
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This result may not be optimal. We expect that under the conditions of Theorem 14.21 one 
can at show that N(E) < C a \E — E \ a near E Q for arbitrary a > 0. But, as long as the 
distribution /i is chosen symmetric and not too small at ±d max , one does not get Lifshits tail 
decay as in ( 1271) . To make this precise, define the Lifshits exponent 7 at Eq by 

= M-log^)) (41) 
' EIE \og(E-E ) v ; 

whenever this limit exists. Note that 7 < 0. If 7 < 0, then it determines the asymptotics of 
the IDS in the sense that, up to logarithmic corrections, N(E) ~ Ciexp(— Ci{E — -Eo) 7 ) as 

E I E . 

Theorem 4.3. Assume that the distribution \x is symmetric and satisfies 

f^i\d max 6, d max j U [ d max , d max -\- c]) ^ C\t (^2) 

for some positive C\ and N and all e > 0. Also assume that the single-site potential q is 
uniformly holder continuous, i.e. that \q(x) — q{y)\ < C^z — y\ p for some C2 and p > and 
all x, y. 

Then 7 = 0. 

In the following we sketch the proofs of Theorems 14. 21 and 14. 31 referring for additional details 
to ffl. 

To prove Theorem \4-'^\ we follow the general strategy of the proof of Theorem 14.11 starting 
with (13Ti) and using the test function ijj u := 0lUn/\\9l u n\\- However, the construction of 
needs to be modified as follows: On each interval [i — 1/2, i + 1/2], i e {1, . . . , L}, un is chosen 
to coincide with a constant multiple of the positive ground state of the Neumann problem for 
— d 2 / dx 2 + q(x — i — uji) on [i — 1/2, i+1/2]. Scaling constants are chosen such that mat (1/2) = 1 
and mat is continuously differentiable throughout [1/2, L+ 1/2]. As we now generally have 
E (uJi) 7^ £xu this leads to extra terms in the bound 

n(e) > ^(\^L^-e'M\ + 2\{e L ^e'M\ 



l \ \\e L ^ 



+ £(Eo(u>i) ~ E ) l ^f 2 <E-E ), (43) 



here := JtiM- 

Due to the symmetry of fi, the numbers log un (i + 1/2), i = 1, . . . , L, are still a symmetric 
random walk (but not Bernoulli) . Versions of the reflection principle and central limit theorem 
for general symmetric random walks and a choice of L as in (13T?]) (with a suitable positive 
constant replacing logr) lead to the bound 



N(E) > ( > jEoiui) - E ) Ij^jj; < E - E 




\0lA 
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> -P > EnUi )-En )<E-En 



> ^ ( ¥(E (ui) - E < ^— ^) ) . ill) 



12 BAKER, LOSS, AND STOLZ 

In Lemma 2.1 of [2 J the continuity of E (-) was shown. The proof given there provides the 
bound 

\E (ai) - E (a 2 )\ p < C J \q(x - a x ) - q(x - a 2 )\ p dx 

for any p > 2. Uniform holder continuity of q gives holder continuity of Eq(-). Using E = 
E (d max ) = E (-d max ) and (H"*} we see that P(£ M - < 5) > C 1 (5/C) N / f) . We plug into 

From this bound, having chosen L through ( |39l) . a calculation shows that the Lifshits exponent 
vanishes. 

The proof of Theorem 14.21 is based on the standard upper bound, e.g. [3], 

N(E)<CF(E Q (H» L )<E). (45) 

We choose L through 

Soe _ 2Cl( L + i) < E — E < So e- 2Cli (46) 
with constants s and G\ to be determined later. By the calculation done in (J7j), 
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Jl/2 l^-l 2 



1=1 



where is the ground state of By a priori bounds (e.g. (16]) there exists C\ > such 

that 

e" Cli < / l^dx < e ClL 
Ji-i 

uniformly in L 6 N, % 6 {1, . . . , L} and all configurations u. Using this C\ in (j46l we further 
estimate 



P(^ « L ) < E) < P ^ < So ) 



< e 



-70 L 



Here the last step is a large deviations bound, which is applicable with suitably chosen s > 
and 70 > due to the assumption (14*01) . Note for this that E(ui) — E are non-negative random 
variables which are strictly positive with positive probability. With this So in (1461) if follows 
that e _7oL < C(E — E ) a , where a := sq/AC\. This completes the proof. 



5. Concluding remarks 

With the above results we have only started to touch the various possibilities for the low- 
energy asymptotics of the IDS in the random displacement model. There are several other 
regimes which we haven't considered yet: 

(i) For one-dimensional random displacement models with non- symmetric distribution, in 
particular the case \i = pSd max + (1 — p)S-d max with p ^ 1/2 we expect that the IDS might 
have Lifshits tails. 

(ii) It would be most interesting to decide if the uniqueness of the minimizing periodic 
configuration established in Theorem 12.4( b) leads to Lifshits tails of the IDS at Eq for the 
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multi-dimensional random displacement model with general (or suitable) distributions \i. Be- 
yond uniqueness of the minimizing configuration this would require to have quantitative results 
on the probability that other configurations have ground state energy near Eq. 

In this context we mention the recent work of Klopp and Nakamura [H] on sign-indefinite 
Anderson models, where some phenomena similar to those found by us for the random dis- 
placement model appear. In particular, they find that Lifshits tail as well as van Hove asymp- 
totics of the IDS at the bottom of the spectrum are both possible in their model, depending 
on the choice of single-site potential and distribution of the random parameters. They have 
informed us about work in preparation [12] which, when combined with the uniqueness re- 
sult Theorem 12.4( b) above, should indeed lead to Lifshits-type asymptotics of the IDS for 
multi-dimensional random displacement models as considered here. This will need suitable 
assumptions on the distribution \x of the displacements, namely that fi is concentrated on the 
corners of [-d max , d max ] d . 

(iii) We have used the non-overlap condition supp \i C [—d max ,d max ] d , d max + r = 1/2, 
mostly for technical reasons. In particular, it is crucial for the Neumann-bracketing argu- 
ments used in [2] and also Section [3] above. However, relaxing this condition will also lead to 
different phenomena. We mention the recent work by Fukushima [8] who studies the random 
displacement model ([I]), (J2J) for positive q and displacements with unbounded distribution 
/i. In this case it is easily seen that the almost sure spectrum is [0, oo), due to the presence 
of large empty regions in typical single-site configurations (while in our setting the spectral 
minimum would be strictly positive). Under this condition Fukushima establishes Lifshits 
tails of the IDS at 0. Another interesting task would be to look at intermediate cases, where 
supp \i is bounded but not small, allowing overlapping finite clusters of single-site potentials, 
but no large empty regions. 

(iv) Under alternative (ii) all random configurations give the same ground state energy. 
This is an example of a random operator with a stable spectral boundary (as opposed to 
fluctuation boundaries). In other examples of this type, for a discussion see Sections 6B and 9 
of [13], this has been found to lead to van Hove behavior of the IDS, i.e. N(E) ~ (E — Eo) d / 2 
as for the unperturbed Laplacian. We also expect this here. 

Acknowledgements: M. L. would like to acknowledge partial support through NSF grant 
DMS-0600037. G. S. was partially supported through NSF grand DMS-0653374. 
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